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Abstract. Consider a probability density in the phase space R.^ X R.^, car- 
ried by the graph of a vector field U m on R", i.e. a Radon measure of the form 
pin = p in (x)5(£, - U in {x)). Let <J> t be a Hamiltonian flow on R N X R N . In 
this paper, we study the structure of the transported measure p(t) := &t#H m 
and of its integral in the § variable denoted p(t). In particular, we (a) provide 
estimates on the number of folds in supp(/i(t)) = <J>t( graph of U m ), (b) estab- 
lish a decomposition of fi(t) into a "regular" component whose integral in the 
£ variable is absolutely continuous with respect to the Lebesgue measure ££ , 
(c) discuss the possibility of atoms for the measure p(t) and (d) construct an 
example in which p(t) is singular with respect to Jf N and diffuse. We conclude 
our study by explaining how our results can be applied to the classical limit 
of the Schrodinger equation by using the formalism of Wigner measures. Our 
results hold for initial momentum profiles If" less regular than C 1 , for which 
the usual notion of caustic is not relevant. The proofs of these results is based 
on the area formula of geometric measure theory. 



1. Introduction 

The subject of this article is the propagation of a certain class of positive Radon 
measures by Hamiltonian flows. 

Let H = H(x,£) be a Hamiltonian of class C 2 on T*R N = Rf x R^. Assume 
that the system of Hamilton's equations 



(1) 



(X t = V s H(X t ,E t ), X (x,0=x, 
\Zt =-V x H(X t ,Et), Ho(z,0=£, 
generates a global flow 

(2) $t(x,£) = (X t (x,£),E t (x,t)). 

Let fi m be a monokinetic measure on R^ x R^ i.e. a Radon measure of the form 

(3) H in &0 = P in (x)6u^ x m 
where U m is a vector field on H N and p m E L 1 (R JV ). Define 

(4) m - **#M m 
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the push-forward of /i in under $ t . Equivalently, fi(t) € C&(R + ; w - M(R N x R w )) 
is the unique weak solution of the Liouville equation 



dtii + {H,n} = 0, 



(5) 

"\t=o 

where {•, •} designates the Poisson bracket 

Our purpose is to study the structure of /i(t) and to deduce from it some information 
on its support 

supp( M (*)) C A t := <f t ({(y, U in (y)) \ y e R"}) . 

While At is the image under $t of the graph of U m , it is not in general the graph 
of a vector field on R N for all values of t. 

When U m — VS m is a smooth gradient field on R", then A t is the union of 
graphs of x i-> V x Sj(t,x), where Sj is a solution of the Hamilton- Jacobi equation 

(6) dtS(t,x)+H(x,V x S(t,x)) = 

defined on some open set of Rt x H N . The graphs of V x Sj are glued along sub- 
manifolds of A 4 where the restriction of the canonical projection 

II : T*R N = R^ x Rf 3 (x, £) ^ x e R W 

is not smooth: see §8 in [5] and §46-47 in [6]. A natural question is to compute, 
or at least estimate, the number of solutions Sj of the Hamilton- Jacobi equation 
needed to obtain A t . 

Equivalently, the restriction to A t of the canonical projection n is in general not 
one-to-one for all t. The question above reduces to the following problem. 

Problem A. To estimate the number Af(t, x) of points in A t whose image under 
n is x. 

The function Af describes the number of folds in A t induced by the Hamiltonian 
dynamics $ t on the graph of U m — even when U ln is not a gradient field. 
Another natural question is the following one. 

Problem B. To study the structure of /i(t) and of its push-forward p(t) under the 
canonical projection n — equivalently, of the first marginal of the measure /i(t) in 
the product space R^ x R^. 

As we shall see below, both problems are intimately related. 

The mathematical problem described above appears in a great variety of con- 
texts. It appeared first in the theory of geometric optics, in the works of Fermat 
and Huygcns: see for instance chapter VII in [15], chapter III in [11] and chapter 12 
§2 in [13]. It appears in the classical limit of quantum mechanics: see for instance 
chapter VII in [16] or [25]. 

Specifically, we shall discuss the following problem in connection with Problems 
A-B: 

Problem C. To study the propagation under the dynamics of the Schrodingcr 
equation of a WKB type ansatz, i.e. of a wave function of the form 

ip in {x) =a ln {x)e lSm(x)/h 
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in the classical limit, in the case of rough (non C 2 ) phase functions S m and square 
integrable (not necessarily continuous) amplitudes a m . 

We recall that the usual approach to the classical limit of the Schrodinger equa- 
tion involves the construction of a Fourier integral type parametrix and applying 
the stationary phase method: see for instance [17] and the proof of Proposition 7.1 
below in the appendix. This obviously requires stringent regularity conditions on 
both the initial phase and amplitude functions in the WKB ansatz. 

In both examples above, the fact that U m is a gradient field is important. There 
are however other types of physical models leading to the same mathematical prob- 
lem even when U m is not a gradient field. Indeed, the Maxwell distribution with 
density p, bulk velocity U and temperature 9, i.e. 

P( x ) r -\$-U(x)\ 2 /26 

converges weakly to the monokinetic measure p(x)6u( x )(€) as # ^ + . The propa- 
gation of this class of measures by the flow $ t generated by the free Hamiltonian 
H(x,£) := ^|£| 2 , i.e. $ t (x, £) = (x + <£,£) can be viewed as the kinetic theory of 
pressureless gases, and appears for instance in a cosmological model due to Zel- 
dovich [27, 9, 8, 7]. When U m is a gradient field, the Liouville equation (5) can 
therefore be viewed as a kinetic formulation of the Hamilton- Jacobi equation (6). 

While the classical mathematical theory of geometric optics or of the semi- 
classical limit of quantum mechanics is centered on the geometry of At in the case 
where both the Hamiltonian H and U m = VS m are smooth, our approach of the 
mathematical problem stated above is centered on the propagation of the monoki- 
netic measure p(x)5jjin r x \ (£) by the Hamiltonian flow. Besides, our analysis on the 
propagation problem is focussed on mathematical methods and results in which 
the initial momentum profile U m is not required to be everywhere differentiable 
— so that A t is not even a C 1 -manifold. Likewise, the Hamiltonian nature of the 
dynamics is of limited importance in our analysis, and there is no need for the mo- 
mentum profile U m to be a gradient — or, equivalently, for A t to be a Lagrangian 
submanifold of T*H N . 

However, we assume that the vector field generating the dynamics is at least 
of class C 1 — or equivalently that the Hamiltonian H is at least of class C 2 on 
x — so that the existence, uniqueness and regularity of the flow $t results 
from the classical Cauchy-Lipschitz theory. The classical limit of quantum dynam- 
ics with rough potentials, for which the existence and uniqueness of the classical 
Hamiltonian flow does not follow from the Cauchy-Lipschitz theory, has been re- 
cently studied in [2]. Our viewpoint in the present paper is different and in some 
sense complementary: we focus our attention to the special class of monokinetic 
measures and to their propagation by smooth Hamiltonian flows, but obtain de- 
tailed information on the structure of the propagated measure p,{t) even for rough 
initial momentum profiles U m . 

The outline of the paper is as follows. In section 2 are collected our assumptions 
on the Hamiltonian H with some elementary estimates on the flow $ t that are 
crucial in the sequel. Section 3 is focussed on the problem of estimating the number 
of folds in the support At of the propagated measure p,(t); our answers to Problem 
A are to be found in Proposition 3.3 d)-e) and in Theorem 3.4 c). The structure 
of fj,(t) itself is studied in section 5; see Theorem 5.1 for our answers to Problem 
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B. Section 4 gathers together a few examples showing that the results in section 3 
are sharp. In section 6, we study different exceptional sets that appear naturally 
in connection with the structure of the projected measure p(t) := H#/j,(t), and 
explain how these sets are related to the traditional notion of "caustic" (introduced 
by Tschirnhaus in the 17th century in the context of geometric optics). Section 7 
discusses applications of the theory presented in sections 3-5, to the classical limit 
of quantum mechanics. While we have not completely solved Problem C, partial 
answers to this problem are stated as Proposition 7.4 b) and c). 



2. On the Hamiltonian flow 

Let H = H(x,£) e R be a C 2 function on R N x R w satisfying the following 
assumptions: there exists n > 0, and a function h e C(R;R + ) that is sublinear at 
infinity, i.e. 

-> as r -> +oo 

r 

such that 

\V ( H(x,0\ <«(1 + |€|) 
(7) \V x H(x,Q\<h(\x\) + K\Z\ 

\V 2 H(x,t)\< K 

for all (x,£) e R w x H N . 

Lemma 2.1. Any Hamiltonian H G C 2 (R W x R w ) satisfying (7) generates a global 
Hamiltonian flow $ t on R N x R^ . The map 

R x R w x R w 9 (t, x, i ^ $ t (.x, e R w x R w 

is o/ dass C 1 . Moreover, for each rj > 0, there exists C n > such that 

sup |X t (a;,0-z| <C„(l + |£|)+»7|a;| 

\t\<T 



for each x, £ E H N , 

\DQ t (x,Z) - Id R « xR « | <e K l 4 l-l 

for all t e R. 

Proof. By (7), one has the following a priori estimates, with the notation (2). 
First 

\X t (x,0 -x\ < f \V 6 H($ s (x,t))\ds < Kt + K [ \E s (x,0\ds 
Jo Jo 



|£ t (*,OI <l£l + I \V x H(<S> s (x,Z))\ds 
Jo 



la 

ft 



and 



<|£|+k J \E s ( Xl O\ds+ I h(X a (x,Q)da. 
By Gronwall's inequality, for all < s < t 

(8) is.fooi < + M*t(*,4)Wt ) ■• 



MONOKINETIC MEASURES WITH ROUGH MOMENTUM PROFILES 



5 



SO that 



\X t (x,0~x\<Kt + K^ eKSds (\Z\+ J Q h(X T (x,£))dT 
< Kt + e Kt (\Z\+f h(X T (x,0)d 



Since h is sublinear at infinity, we have, for every R > 

h(r) 

(9) h(r) < l [0 ,ffl(r) sup h(r) + l( R)+00 ) (r)r sup r)<M R + rm R , 

0<r<R r>R r 



where 



so that 



Therefore 



h(r) 

Mr = sup h(r) and m R = sup 

0<r<fl r>R r 



m R — >• as R — > +oo . 



Jo 

< (K + M K e Kt )t + e Kt |e| +m K e Kt |a;| + m fl e Kt f \X a (x,£) - x\da . 

Jo 

By Gronwall's inequality, 
(10) 

\X t (x,0 -x\< ((« + M fl e Kt )t + e Kt |e| + ro^'larDe*"* 6 "* 

< Kte*" 1 " 6 "* +M K te t(K+m « eKt) + \£\e t( - K+mReKt '> + m R \x\e t{K+mReKt) . 

The same estimates hold for — T < t < after substituting \t\ to t. 

In view of (10)-(8), for each (x,£) £ R N x R N , the trajectory (x,£) ^ $t(x,€) 
cannot escape to infinity in finite time, and is therefore globally defined. 

Besides, since m R — > as R —¥ +oo, the estimate (10) obviously implies the first 
inequality in the lemma with 

7] := m R e T ^ +mRe ^ and C„ := (1 + kT + M R T)e T ^ +mRe " T ^ . 

Since H e C 2 (R N x R w ), the map {t,x,Q ^ $t(x,0 is of class C 1 on its 
domain of definition R x R N x R N . Differentiating the Hamilton equations (1) 
with respect to the initial condition, one finds that 

f D'X t = +V^tf ($ t ) ■ DX t + V^iJ(* t ) ■ DE t , 
{ DE t = -VS, x H(*t) • DX t - Vl ( H(4> t ) ■ DE t , 

so that 

r\t\ 

\D$ t - IdRNxRN | < k / \D$ s \ds. 

Jo 

The second inequality in the lemma follows from Gronwall's inequality. □ 
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3. On the number of folds in A t 

Let U in e C(K N ;R N ) satisfy the condition 

\U m (v)\ 

(11) 1 , \ gn -» asM->0. 

1 3/1 

The present section is focussed on the structure of the set A t in the introduction, 
defined as 

(12) A t :=M{(y,U in (y))\yeR N }), teR. 
Consider the map 

(13) F t : R N 3 y ^ F t (y) = X t (y, U m (y)) e R N . 

Whenever U m is differentiable at y, the map F t is also differentiable at y by the 
chain rule; for any such y, define the absolute value of the Jacobian determinant 

(14) J t (y) = \det(DF t (y))\. 
We also introduce the set 

C := {(t,x) e R x n N I F t -\{ X }) n ^({O}) + 0} , 

(15) „ 
C t :={xeR N \(t,x)eC}. 



For lack of a better terminology and by analogy with geometric optics, C will be 
referred to as the "caustic" set. 

Proposition 3.1. Under the conditions above, for each (eR, the map F t is proper 
and onto. Moreover 

fta\ \ F t(y)-y\ „ , , 

(16) sup 1— j > as \y\ -> oo . 

\t\<T \y\ 

The proof of this proposition and the next one will use the following topological 
argument. 

Lemma 3.2. Let g : R N — > R N be a continuous map satisfying the following 
condition: for some R > 

(g(x)\x) > for all x e R w such that \x\ = R. 

Then 

a) there exists x G R w such that \x\ < R and g(x) = 0; 

b) if g is of class C 1 on R N and is a regular value of g, then 5 _1 ({0}) n B(0, R) 
is finite and #(.g _1 ({0}) n B(0, R)) is odd. 

Proof. Consider the homotopy G e C([0, 1] x R w ; K N ) defined by 

G(t,x) = tx + (1 - t)g(x) . 

One has 

G(t, x) ^ whenever t e [0, 1] and \x\= R. 

Indeed, G(l,x) = x ^ if |a;| = R > 0; besides, if t e [0, 1[ and G(t,x) = 0, one 
has 

g(x) = ~Y—j. x 80 that i.9{x)\x) = - — -\x\ 2 = -Y— t R2 < 
for all x e R N such that |x| = R, which contradicts our assumption. 
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By the homotopy invariance of the degree (see Properties 7, 8 and Theorem 12.7 
in chapter 12, § A of [24]) 

d(g,B(0,R),0) = d(I,B(0,R),0) = l. 

This implies a). Moreover, if g is of class C 1 on H N and is a regular value of g, 
all the elements of g -1 ({0}) are isolated points by the implicit function theorem, 
so ff _1 ({0}) n B(0,R) is finite. Besides (see Property 2 in chapter 12, §A of [24]) 

d(g, B(0, R), 0) = J2 sign(dct Dg{x)) = 1 . 

xeg- 1 ({0})r\B(G,R) 

Therefore, there exists an integer m G N such that 

#{x e B(0, R) | g(x) = and det(Dg(x)) > 0} = m + 1 
#{x e S(0, i?) | 5(x) = and det(Dg(x)) < 0} = m 

so that #( 9 - 1 ({0}) n B(0,R)) = 2m + 1, which proves b). □ 

Proof of Proposition 3.1. By the first inequality in Lemma 2.1 and the condition 
(11 on ?7 m , for each r\ > 0, one has 

— |F t (y)-y| 

hm sup — : < r] , 

ii/K+oo |t|<T \y\ 

so that (16) holds. 

Because of (16), the continuous map F t satisfies 

(17) (F t (y)-x\y) = \y\ 2 + o(\y\ 2 ) as |y| ^ +oo 

so that F t is onto by applying Lemma 3.2 to the map g : y M> F t (y) — x. On the 
other hand 

\F t {y)\ -> +oo as \y\ -> +oo 
so that is proper. □ 

By Proposition 3.1, we know that, for each (t, x) e R x R w , the equation 
F t (y) = x has at least one solution y G when U ln is a continuous vector field 
sublinear at infinity, i.e. satisfying (11). In the next proposition, we study the 
number Af(t, x) of solutions of this equation in the case where U m is of class C 1 at 
least. Equivalently, Af(t,x) is the number of intersections of the manifold A t with 
T*H N ~ {x} x R w . Therefore, the integer-valued function Af measures the number 
of folds in the manifold A t resulting from the interaction of the Hamiltonian flow 
$ t with the initial profile U ln . 

Proposition 3.3. [Smooth case] Assume that (11) holds for U in G C^R^.R^) 
and that 

\DU m (y)\=0(l) as\y\^ +X1 . 

a) For each t € R, one has £? N (C t ) = 0. 

b) The set C is closed in R x R w . 

c) For each (t, x) e R x R N \ C, the set Ff ({x}) is finite, and henceforth denoted 
by 

{Vj(t,x), j = l,...,Af(t,x)} . 
The integer N is a constant function of (t, x) in each connected component of 
R x H N \ C and, for each j > 1, the map yj is of class C 1 on each connected 
component of R x R N \ C where N > j. 
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d) There exists a < < b such that C(~l((a, b) x K N ) = and M = 1 on (a, b) x R N . 

e) Af(t, x) is odd for each (t, x) e R x H N \ C. 

Proof. If U in e C^R"), the map F t is of class C7 1 from R w to R w , being the com- 
position of C 1 maps. Since Ct is the set of critical values of Ft, one has ££ N {Ct) = 
by Sard's Theorem (in the equal dimension case), which proves a). 

Pick any sequence (t n , x n ) € C such that (t n ,x n ) — > (t, x) in R x R N as n — > oo. 
By definition of C, there exists y n such that Ft n (y n ) = x n . Since the sequences 
i„ and x„ converge, they are both bounded. Let T = sup„ \t n \; assume that the 
sequence y n is unbounded; if so there exists a subsequence y ni such that \y ni | — > oo. 
By (16) 

|a;„, = \F t (y ni ) - y m \ < sup \F t (y ni ) -y n ,\ = o(\y„,\) 

\t\<T 

so that \x n , | ~ \y ni \ as n; — > +oo. Since this contradicts the fact that the sequence 
a; n is bounded, we conclude that the sequence y„ is bounded. Therefore, there 
exists a convergent subsequence y nk of y„; call y its limit as — > +oo. Passing to 
the limit in both relations 

*t»„(»nj = ^n fc and J tnk (y nk ) = |det(VF tiifc (y„J)| =0, 

we conclude that (t, x, y) satisfies 

F t (y) = x and J t (y) = 

and therefore that (t, x) G C , which proves b). 

For (t, x) e R x R^ \ C, all the solutions of the equation F t (y) — x = are 
isolated by the implicit function theorem. The set of all such points, F t _1 ({a;}), is 
therefore compact since F t is proper. By the implicit function theorem, the integer 
AT is a locally constant function of (t, x) <E R x H N \ C, and is therefore constant on 
each connected component of R x R^ \ C. Let j <G N*, and let O be a connected 
component of R x R N \ C; by the implicit function theorem yj € C 1 (il). This 
proves c). 

Assume mf{t > | C t ^ 0} = 0. Then, there exists (t n ,x n ,y n ) such that 

t n -> + , F tn (y n ) = x n , and J tn (y n ) = . 

Assume that a subsequence y nk of the sequence y n is bounded. Up to extraction 
of a subsequence, one can assume that y nk -4- y, so that = Jt nk (yn k ) ~> Jo(y)- 
But since F t — Id R jv, one has Jo(y) = 1. Therefore \y n \ — > +oo. By the second 
inequality in Lemma 2.1 

\D x X tn (y n , U m (y n )) - Id R « | < e K ^ - 1 , 
\D^X tn (y n , U m (y n )) ■ DU m (y n )\ = O (e^\ - l) , 

so that 

= JtM = I det(D x X tn (y n , U m (y n )) + D^ B (y n , U m (y n )) ■ DU m (y n ))\ 
— > | dct(Id R «)| = 1 as n — > oo . 
Thus the assumption i„ — > leads to a contradiction. Therefore, 

inf{i >O|C7 t ^0} = 6>O. 

By the same token, 

sup{i <O|C t ^0} = a<O. 
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Thus (a, 6)xR JV is contained in the connected component of {OjxR^ in RxR JV \C. 
Since Fq = Id R iv, one has Af(0,x) = 1 for all x € R w , and since Af is constant on 
each connected component of R x R w \ C, one concludes TV = 1 on (a, b) x R w , 
which proves d). 

If (t, x) £ R x R N \ C, the point a; is a regular value of F t . Since F t is proper, 
Ff ({x}) is compact, and therefore bounded. Pick R > such that 

F^dx^cB^R). 

By (17) and Lemma 3.2 b) applied to the map g : y F t (y) — x 

#F t -\{x}) = #(F t - 1 (M) n B(0,R)) is odd, 

which proves e). □ 

When U m is not of class C 1 , the arguments used to prove Proposition 3.3 are no 
longer valid. However one can still obtain some information on the number Af(t, x) 
of solutions y of the equation F t (y) — x by a completely different method, involving 
the area — or co-area formula. 

Assume that U m G C(R N ; H N ) satisfies (11) and that its gradient (in the sense 
of distributions) DU ln satisfies the condition 

(18) d l U k n \ Q e L N '\{1) for each bounded open n C R w , 

for all k, I = 1, . . . , N. We recall that a measurable function / : ft — > R belongs to 
the Lorentz space L Ar,1 (il) if 



poo 

/ (j? N ({xen\\f(x)\ > \})) 1/N d\< 
Jo 



oo . 



By Theorem B in [14], the vector field U m is diffcrentiable a.c. on H N . Let E be 
the ££ ^-negligible set defined as 

(19) E := {y £ K N \ U m is not diffcrentiable at y} . 

By the chain rule, the absolute value of the Jacobian determinant 

My) = I det(D x X t (y, U m (y)) + D s X t (y, U in (y))DU in (y))\ 
is defined for all (t, y) £ R x {K N \ E) . 

Henceforth, the notation J t _1 ({0}) designates the set 

(21) J^({0}):={y£R N \E\J t (y)=0}, 
and we shall also consider the set 

(22) Z t :=J-\{0})UE. 



Theorem 3.4. [Rough case] Assume that U m £ C(R W ;R W ) satisfies (11) and 
(18). 

a) for all t £ R, the function J t £ Lj oc (R N ); 

b) for all tgR, there are finitely many solutions y of the equation 
(23) F t (y) = x 

for a.e. x £ H N — in other words, Af(t, x) is finite for a.e. x £ H N , for all t £ R; 



10 C. BARDOS, F. GOLSE, P. MARKOWICH, AND T. PAUL 

c) for each bounded B C R w , denote by x) the number of solutions y € B of 

(23) ; then, whenever B is bounded, for each teR and each ueN, one has 

% N ({x e R, N \Af B (t,x) >n})<- f J t (y)dy 

n J B 

< ±e N ^\\l + (l-e- K W)\\DU in \\% N{B y, 

d) let t e R; i/ien /or a.e. x <G R^, aZZ solutions y of (23) satisfy 

yeR N \E and J t (y) >0; 

ej /or allT>0 

Jf\{(t, y) G [— T, T] x R N | F t (y) = x}) < +ao 

/or a.e. i£ R w . 
Remarks. 

a) By the first statement in Proposition 3.1, J\f(t,x) — A/rat (i,x) > 1 for all 
(t,x) e Rx R N . 

b) Even in the smooth case, i.e. assuming in addition that U m € C7 1 (R Ar ), state- 
ment c) in Theorem 3.4 provides information on the number of folds of A t , that is 
the image under $ t of the graph of U m , which seems to be new at the time of this 
writing. 

Proof. By the second estimate in Lemma 2.1 

(24) \D x X t (y, U m (y)) + D t X t (y, U m (y)) ■ DU m (y)\ < e K ^ + (e K ^ - l)\DU m (y)\ 
so that 

Mv) = I det(D x X t (y, U m (y)) + D c X t {y, U m (y)) ■ DU m (y))\ 
< e JV«l*l(i + {i- e -«-\t\)\DU in {y)\) N 

by Hadamard's inequality. Since U m satishes (18) and since L N ' l {0) C L N (fl), 
this inequality implies statement a). 

Since the map F t is proper by Proposition 3.1, the set K ti R = F t " 1 (B(0, R)) is 
compact for each R > 0. By the area formula (see Theorem 3.4 in [20] and Theorem 
A in [14]) 



L 



#(F t - l {{x})C\K t , R )dx= [ J t (y)dy < +oo . 

■Ik, n 



R™ JK t , R 

Therefore #F t _1 ({x}) < oo for a.e. x e B(0,R); since this is true for all R e N, 
one concludes that f^F 1 T l ({x}) < oo for a.e. x e R w , so that b) holds. 

Let B be a measurable subset of R w ; applying again the area formula shows 
that 

/ N B {t,x)dx= / J t {y)dy. 
By the Bienayme-Chebyshev inequality, for each n > 1 

^({ieR^|^(t,i)>fi})<-/ My)dy, 

n Jb 

which is precisely the first inequality in c). The second inequality follows from (25) 
and Holder's inequality. 
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Let n = 1 and Z t := J t ~ ({0}) U £7. Let £> = Z 4 n i^t,i?; the set £? is measurable 
and bounded. Then 

/ J t {y)dy = 
Jb 

since J t (y) = for all y e B \ E and £? N (E) = 0. Applying c) shows that 

& N ({x eR N \ \x\ < R and M Zt (t, x) > 1}) = , 

which entails d) by monotone convergence, letting i? € N tend to infinity. 
Consider next the continuous map 

F: [-T,T]xK N 3 (t,y) ^ F(t,y) e K N . 

In view of (16), \F(t, y)\ — > oo as \y\ — > +oo uniformly in f e [— T,T]. Therefore, 
the set := ^"^(O, R)) is compact for each R > 0. Then, for each f e [-T, T] 
and each y € R w \ E, the Jacobian DF(t, y) is the column- wise partitioned matrix 

DF(t,y) = [V(t,y),M(t,y)], 

with 

F(t,y) = V c £f(<i> t (y,[/ m (y))) 

and 

M(t, y) := D x X t (y, t/ m (y)) + £> £ X t (y, XP^DlP^y) . 

Therefore, 

DF(t, y)DF(t, y) T = K(t, y)F(t, y) T + M(i, y)M(t, y) T 
so that, by the co-area formula (Theorem 1.3 in [21]) 

/ JP x {F- x {{x})nK R )dx 

JR N 

= f J det( V(t, y)V(t, y) T + M(t, y)M(t, y) T )dtdy . 
Jk r v 

By Lemma 2.1, (t, x, £) i-> <&t(:r, £) is of class C 1 on R x R N x R w , so that the map 
(t, y) i-> V(i, y) is continuous on R x R w , and therefore bounded on the compact 
Kr. On the other hand, by (24) 

sup \M(t, y)\ < e KT + (e« T - l)|DIT n (i/)| € Lj^R") , 

|t|<T 

since t/ m satisfies (18). Denoting 

:= {y e R w | there exists t e [-T, T] s.t. (i, y) e K R } 

that is compact in R N (being the projection of the compact Kr on the second 
factor in R x R w ), one has 

\\VV T + MM t \\%% Kr) < 2 N ^\\V\\^ {Ka ^ N +\K R ) + 2^T\\M\\ n lN[KIr) 

< oo . 

Therefore Jrf? 1 (F^ 1 ({x})<~)K r) < +oo is finite for a.e. x <G B(0, R), and since this 
is true for all R <G N, one concludes that 1 (F~ 1 ({x})) < +oo for a.e. x G R ff , 
which is statement e). □ 
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Theorem 3.4 suggests considering the sets 
(26) C' t = {x E~R N \ F^dx}) is infinite} 

for all t e R. 

On the other hand, the definition of the caustic (15) in the smooth case should 
be modified as follows when the continuous vector field U ln is not of class C 1 but 
satisfies (11) and (18): 

C t :={xeR N \F t - 1 ({x})riZ t ^0}, 

W C:={)C t , 

ten 

where Z t is defined in (22). Notice that this definition coincides with (15) whenever 
U m is differentiable everywhere. However we do not know whether C or C t are 
closed whenever U m is not of class C 1 . 

The first part of statement c) in Proposition 3.3 is equivalent to the inclusion 

(28) C' t C C t for all f e R . 
Statements b) and d) in Theorem 3.4 can be recast as follows 

(29) & N (C' t ) = and Sf N (C t ) = , for all t E R . 

This is the analogue of statement a) in Proposition 3.3 in the case of a rough 
U m . Notice that (29) is a consequence of the area formula while statement a) in 
Proposition 3.3 follows from Sard's theorem (see Remark 2.97 on pp. 103-104 in 
[3], discussing the relation between Sard's theorem and the co-area formula). 

4. Some examples 

While the notion of caustic C naturally occurs in studying the geometry of A t in 
the smooth case, its relation to the number of solutions y of the equation F t (y) = x 
is slightly less obvious, as shown by the following examples. In all these examples, 
the Hamiltonian is 

generating the free flow 

$ t : RxR3(i,()i->(i + ^)eRxR. 

Thus 

F t : R3y^y + tU in (y) eR. 

The first example below shows that the set Ff ({x}) may be finite even if x e Ct- 
In other words, it may happen that the inclusion (28) is strict for all t£R. 
Example 1. Set N = 1, and let U m be real analytic on R and satisfy (11). 
Therefore F t is real analytic on R for each t € R. For each x <G R, the set 
Ff ({x}) is the set of zeros of the analytic function y i-> F t (y) — x. Since U m 
satisfies (11), one has F t (y) <~ y as \y\ — > oo, so that, for each fgR, the function 
F t is not a constant. Therefore the zeros of y F t (y) — x are isolated for each 
(i, x) € R x R. Equivalently, the set F^ 1 {{x}) consists of isolated points. On the 
other hand, for each (t,x) € R 2 , the set F t " 1 ({x}) is compact since F t is proper by 
Proposition 3.1. Therefore the set F t _1 ({x}) is finite for all (t, x) € R 2 . 

The next example shows that the set F 1 T l ({x}) may be infinite for infinitely 
many times t — in this case, for all t in a nonempty, open interval of R. 
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Example 2. Set N = 1, and let U m be denned by 

{tanh(z) sin(ln \z\) if z > , 
ifz = 0. 

Clearly U in G C X (R*) and, for all z G R*, one has 

(C/ m )'(z) = (1 - tanh 2 (z)) sin(ln \z\) + tanh ( z ) cos (ln |z| 



Observe that (U m )'(z) = s/2 sin(ln \z\ + § ) + 0(z 2 ) does not have a limit for z -> 0, 
so that C/ m ^ C^R). On the other hand 

sup \U m (z)\ = 1 and sup \(U m )'(z)\ < 2 

so that U m is Lipschitz continuous on R and the map F t is proper on R. 

Let x = 0; for each t such that |i| > 1, the set F 1 T l ({0}) is infinite. 

Indeed, the set Ff 1 ({0}) obviously contains z — 0; besides, F i T 1 ({0}) contains 
also a decreasing sequence y n — > as n — > +oo, which satisfies 

sin(ln \y„\) — > — l/t as n — > oo . 

Therefore F 1 ~ l ({0}) is countably infinite whenever \t\ > 1. 
Notice that, in example 2, whenever T > 1 

{(*, y) G [-T, T] x R | F t (y) = 0} = ([— T, -1] U [1, T]) x ({0} U {y n \ n > 1}) 
so that 

Jf 1 ({(t, y) G [-T, T] x R | F t (y) = 0}) = +^ . 
In other words, the set of points x for which 

^ ({(*, y) G [-r, T] x R | F t (y) - a;}) = +^ , 
while S£ x -negligible by statement e) of Theorem 3.4, may be non empty. 

There is another interesting observation in connection with Example 2. Observe 
that G C[ and that 

^({O}) = {0} U {y n | n > 1} , whenever |i| > 1 . 
Nevertheless, for each n > 1 



(U m )'(Vn) = sin(ln \y n \) + cos(ln \y n \) + 0(\y n \ 2 ) -+ ~ + ^1-1^-1 

if |t| > 1, and F t is not differentiable at y = 0. Thus y = is in the exceptional set 
E where the Lipschitz continuous function U m is not differentiable, so that £ Cj 
with the definition (27) although 

F t - 1 ({o})n.j t - 1 ({o}) = 0, 

so that would not belong to C t had we kept the classical definition (15) in the 
case of non everywhere differentiable U m profiles. 

The next two examples show that F t -1 ({x}) can even be uncountably infinite, 
even for a smooth profile U ln . 
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Example 3. Set N = 1, and let U m be defined by 

' + 1 if z < -1 , 

U m (z) :-- 



— z 

- 1 



if \z\ < 1 . 
if z > 1 . 



Consider the equation 

F t (y) :=y + tU m {y) = x 
with unknown y; for t < 1, its solution is unique and given by 

( x t if x < t — 1 , 



if |x| < l-i, 
if x > 1 - £ . 



For t — 1, the solution is 



2/= S 



For f > 1, the solution is 
' x — t 



y= \ 



x — t . 
x + t 



( x-t a x <o, 

any y £ [-1,1] if x = , 
x + 1 if .x > . 



if x < 1 - t , 
and x + t ifl — t < x < t — 1, 
if x > t — 1 . 



1-i 



Now is finite for all x G R whenever t ^ 1, while Ff ^{O}) = [-1, 1]. 

Example 4. In the previous example, U m is Lipschitz continuous but not of class 
C 1 . Yet the same phenomenon can be observed by smoothing U m near z = ±1. 
Regularize J7 in and obtain C/ e in G C°°(R) so that 

supp(C/™ - U m ) C [-1 - e, -1 + e] U [1 - e, 1 + e] , and 

jjin < jjin Qn [_ 1 _ £j _j + e ] and jjin < on [1 _ e , 1 + e ] . 

In that case 

F 1 - 1 ({0}) = [-l + e,l-e]. 

Indeed, since U™ < U m on [-1 - e, — 1 + e], all the points on the graph of [7 e m with 
abscissa in (—1 — e, — 1 + e) will reach x = after t = 1. The same is true of the 
points on the graph of U l t n with abscissa in (1 — e, 1 + e). Thus the regularization 
does not affect the dynamics of the points with abscissa in (— 1 + e, 1 — e) for all 
t € [0, 1], and in particular for t = 1. 

5. On the structure of /j,(t) 

Throughout this section, we assume that U m G C(R N ; R N ) satisfies the sublin- 
earity condition (11) at infinity and the regularity condition (18). 

Consider a monokinctic measure \x m of the form (3) with p in e L 1 (R JV ), whose 
action on a test function \ •= Cb(R N x R w ) is given by the formula 



>x) == / x(y, 



U in (y))p in (y)dy. 
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In other words 

Let H e C 2 (TL N x K N ) satisfy (7), and let i> t be the Hamiltonian flow generated 
byH. 

For all t € R, let p,(t) = $ t #/i m be the push-forward of p m under <I> t , defined 
as follows: for each test function \ e C&(R W x R w ), 

(30) <ji(t), X ) ■•= (M™ X o $ t ) = / X (* t (y, (7 m (y)))p m (y)rfy . 

Finally, let p(i) be the measure on R w defined as 

(31) p(t) := n#/x(t) 

where II : T*H N ~ R w x R w 9 (ar,£) i-» x e R w is the canonical projection. In 
other words, for each test function <p <= Cb(R N ) 

(32) ( P (t)A)=f d>(X t (y,U in (y))p in (y)dy 
or, equivalently 

(33) P {t) = F t #p in . 
We shall also use the following definition 

(34) P t :={yER N \E\ J t (y) > 0} . 

Our main result in this section bears on the structure of the measures p(t) and 
its projection p(t). 

Theorem 5.1. Under the assumptions above, let p m e L 1 (R N ) be such that p m > 
a.e. on R w . 

a) for each teR, 

p m l Pt ^ N < J t & N , and p m lz t ^ N -L Jt^ N 

b) for each t € R, 

M*) = M*) + A*s(*) /i (i) _L n s {t) , 

w/iere 



[ /i a (i) := $ t #(j2f ® p in (i)lz t W^- (l )) • 
Moreover 

c, ) for each t € R, one /las 

supp(/i(t)) C A t := * t ({(y, (7 m (y)) | y € R^}) ; 

d, ) /or eac/i ieR, ffte measure 

Pa (t) :=n# Ala (t)«^ JV ; 

^(X)= W; (X) £ qj&N) M 

where C' t is defined in (26); 
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e) for each t € R, the measure 

p s (t) := n#/x s (i) is carried by C t ; 

in particular 

p a (t)± P s(t). 

A few remarks are in order before we give the proof of Theorem 5.1. In the 
smooth case — i.e. whenever U m G C 1 (R JV ;R JV ) and satisfies (11), we recall 
that Ct is closed in R N for each t £ R, by statement b) of Proposition 3.3. For 
any given t £ R, let \ = x( x ,£.) be a test function in C^R^ x R N ) such that 
nsupp(x)nC t = 0. By (30) 

(M*),X>= / x(My,U m (y)))p m (y)dy. 

JR™ 

Since II(supp(x)) is compact and included in the open set R w \ Ct, it intersects 
at most finitely many connected components of TL N \ C\ ■ Assume without loss of 
generality that II(supp(x)) is connected, so that it intersect exactly one connected 
component il of R^ \C t ; on 0, the integer- valued function N is a constant de- 
noted by N n , by statement c) of Proposition 3.3. With the notation used in that 
proposition, for each x € II(supp(x)) C O 

F t -\{x}) = { yj (t,x)\j = l,...,M n }, 

and 

</,•(*,•) eC 1 (ft) for all., I V". 

Therefore yj(t, •) is a C 1 -diffcomorphism from f2 on its image Oj, with inverse F t . 
Thus 

F^itt) = (J Oj and n <3 3 ■ = if i 7^ 3 , 
j=i 

so that 

x (t,x)p(t,dxdO = Y, / x(Mv),Zt(v,u in iv)))p in {v)dv- 

f2xR« j = l^°3 

In each of the integrals on the right hand side, F t is a C 1 -diffcomorphism mapping 
Oj on f2, so that, changing variables, we see that 

/ X (F t (y),E t (y,U m (y)))p m (y)dy 
Jo, 

= [ x(x,S t ( % (t,x),(7 m ( % (t,x))))p m ( % (t,x))|det(V ;c% (t,x))| ( ix. 
Jn 

Since 

\det(D x yj(t,x))\ = J t (Vj(t, x))' 1 , 
we conclude that the restriction of p,(t) to (Rx R N )\C is a measure- valued function 
of (t, x) given by the following formula: 

(35) p(t,x, ■):= j t {y^ t ~yj SE t(yj(t,x),u*"( Vj (t,x))) 

whenever (t,x) ^ C. This formula is strikingly similar to the one giving p a {t) in 
statement b) of Theorem 5.1. There are however subtle differences, which we shall 



II 
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discuss in more detail in the next section. At this point, it suffices to say that 
Theorem 5.1 provides a formula for p,(t) that holds globally on x H N instead 
of (R w \ C t ) x R w , and that the argument above requires more regularity on U m 
than assumed in Theorem 5.1. 

Proof of Theorem 5.1. Let A C R w ; then the condition 

J t (y)dy = implies that Jt(y) = for a.e. t/ei. 



I A 

Therefore £g N \P t n A) = so that 



/ (p in l Pt )(y)dy = f p in (y)dy = 0. 
J A Jp t nA 



Thus p in l Pt & N < J t ^ N . 

On the other hand, for each teR 

R w = p t u z t with p t n z t = . 

Since 

J t (y) = for y e Z t \E i.e. ££ N -a.e. on Z t 

while 

p in (y)lz t (»)=0 foraUyePt, 

we conclude that p m lz t ^ N -L J t Jz? N , which proves a). 
Define 

»™:=J? x N ®(p in l Pt )(x)6™ in{x) , 
^■.= J?»®{p in l Zt ){x)5™ n{x) , 

and, for each t e R 

so that one has indeed 

M*) + M.(t) - ^#CC + MD = *t#M m = M*) • 
Then there exists a unique b e L 1 (R 7V ; J t Jz? N ) such that 

by the Radon-Nikodym theorem. Thus, for each \ G ^(R^ x R w ), by the area 
formula (see Theorem 3.4 in [20] and Theorem A in [14]) 

(M«(*),X>= / x(P t (y),S t (y,C/»(y)))6(y)J t (y)dy 

= ^1 E Ky)xfo3t(y,^ n (y)))j<fe 



/ 



In the formula above 



K yeF-\{x}) 

E b (y)(S~ t (y,u^(y)),x(x,-)) } dx. 

K yeF t -\{x}) 



_ d(ri Pt j? N ) 

d(J t J? N ) 
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is the Radon- Nikodym derivative of p ln lp t with respect to Jt^? N ■ Since J t > on 
the set P t 

(36) b = f-—^ a.c. on K N , 

Jt 

which proves b). 

Formula (3) obviously implies that 

suppOOcAo :={(y,U in )\yeR N }. 

Since p(t) = $ t #p m , one has 

supp( M (*)) C * t ({(y, U in ) | y € R w }) - A t , 

which is precisely statement c). 

In view of the first formula in b) and of the definition p a (t) = U#p, a (t) of the 
measure p a (t), one has 



Pa(t) = [ £ b(y) ] if 
yyeF-^lx}) 



N 



with b as in (36). The set F t _1 ({a;}) can obviously be infinite, in which case the 
sum above can be undefined. However, this occurs only if a; € C[ as defined in (26). 
Since C' t is if ^-negligible 

Pa (t) = f t J? N 

with 

f t (x):=l RNXC ,(x)l ]T b(y) 

\l/€F- 1 ({x}) 

Thus, p a {t) is of the form p a {t) = ,ft^ N with f t > measurable on R w . Besides 
< p a (t) < p(t) and / p(t,dx) = / p in {x)dx) < oo 

which implies in particular that p a {t) <C if w with the formula for the Radon- 
Nikodym derivative as in d). 

Consider the measurable set A := H N \ C t - By (33) applied to p m lz t instead of 
p m , one has 



Ps (t)(A)= P m (y)iz t (y)dy= p m (y)dy = 0. 

Jf-\a) JF~\A)nz t 

Indeed, by definition of C' t , one has Ff 1 {A) (~l Z t = 0. In other words, p s (t) is 
carried by C t . 

Finally, since p a (t) <C if N and p s (t) is carried by C t which is if -^-negligible by 
(29), we conclude that p a (t) -L p s {t) which is precisely statement c). □ 
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6. On the caustic and other exceptional sets 



In the case of a smooth U m profile — i.e. when U m G C 1 (R Ar ;R Ar ) satisfies 
(11), the caustic C is the only exceptional set — C t being equivalently defined as 
the image under the projection II of the set of points in the manifold A t in (12) 
where the restriction II I . is not diffcrcntiable. 



When U m is not everywhere diffcrcntiable, this definition of C t does not make 
sense in general since A t is not a differentiable manifold in the first place. In such 
cases, it is more natural to consider the measures p{t) and p{t) instead of the sets 
A t and Ct — all the more so since Ct may not even be closed in R*. Thus, even 
though p a (t) is concentrated on C t , one cannot say that C't is the support of p s (t) as 
C t may not be closed. On the other hand, the inclusion suppp s (i) C C t is of little 
interest as C t might be dense in some domain of R w . Although p s {t) is concentrated 
on Ct, this obviously does not characterize Ct (if a measure is concentrated on a 
set, it is also concentrated on the complement in that set of any negligible set for 
that same measure). 

There are analogous difficulties with the absolutely continuous part of the mea- 
sure p a {t)- In formula (35), the restriction of p a {t) to (R N \C t ) x R N is viewed as 
a function of x G H N \ C t with values in the set of Radon measures in the variable 
£ G R ff . This viewpoint is obviously not appropriate if U m is not at least of class 
C 1 — for instance, if the set Ct is dense in some domain of R w . In Theorem 5.1, 
the measure p is a weakly continuous function of the time variable t with values 
in the space of Radon measures in the variables (x,£), and is therefore globally 
defined on x R^. Obviously, the ratio p m lp t /J t is just one possible choice of 
the Radon-Nikodym derivative d(p m lp t ^ N )/d(Jt^f N ) and could be modified ar- 
bitrarily on any set of Jt Jzf^-measure — which could be of positive J<? ^-measure, 
as in Examples 3-4 above. The difference induced in the expression 



by two different choices of the Radon-Nikodym derivative d(p m lp t N )/d(J t J? N ) 
is of the form 



yeF-\{x}) 

with m > being a Jt-Sf^-measurable function such that m = a.e. on H N \ Z t . 
Therefore 



(37) 







so that 



yeF-'dx}) 



% N - a.e. on R 



and 



yeF-\{x}) 



m(y)S U in {y) = 



as a measure on R x R . 
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On the other hand, even in the smooth case, formula (35) is not enough to 
define completely /j,(t), as it fails to capture singular parts of the measure carried 
by Ct x R w . For instance, in Example 4, for t = 1, if one computes the restriction 
of the measure /z(l) to (R w \ C\) x R N by formula (35), one obtains 

/ (f n{t,x,d0)dx= f f£ P m (y)dy+ f + °° P m (y)dy 



* \ p m (y)dy 



R 



/ (f fi in (x,d0)dx 



unless p m = on [-1 + e, 1 - e] = Jf ^{O}). Since this is obviously incompatible 
with the conservation of the total mass of the measure 

// p(t,dxdO= ff $> t #n in {dxd£)= [[ p m (dxdO 

iVxR" iJR«xR» JJr«xR» 

by the transportation under the flow $ t , this example confirms that (35) alone is 
not enough to compute p(t). 

The expression involved in the formula for p a in Theorem 5.1, i.e. 

- n* n 1 

(38) ]T e—£L iy) s utn(v) 

seF t "'(W) * 

which is analogous to formula (35) a priori makes sense only ii x E ~R N \ C' t . In 
the rough case, i.e. if U m satisfies only the regularity properties in Theorem 3.4 or 
Theorem 5.1, it might on principle happen that 

C' t n (R N \C t )?0. 

Whenever this is the case, the sum in formula (38) above may fail to be defined 
everywhere on R N \C t , as would occur for U m of class C 1 , by statement c) in 
Proposition 3.3, or (28). Since this may occur at most on a Jzf ^-negligible set, the 
expression for p a in Theorem 5.1 is well defined as a measure on R N x H N while 
the one in (38) may not be defined everywhere on H N \ C t as a function of x with 
values in the space of measures in the variable £ <G H N . 

Another strategy is to abandon the idea of viewing the caustic C as an excep- 
tional set unambiguously defined by the measure p(t), and to investigate those 
exceptional sets that can be uniquely defined in terms of this measure. 

This leads for instance to consider the set 

(39) C'l := {x e n N l j? N (F t -\{x}) n J t -\{0}) > 0} , 

which is of particular interest, as shown by the following result. 

Theorem 6.1. Assume that H e C 2 (R Ar x R w ) satisfies the conditions (7), while 
U m G (^(R^jR*) satisfies the sublinearity condition (11) at infinity and the reg- 
ularity condition (18). Then 

a) C'{ C C t n C' t for each t G R; 

b) for each t € R and each x <G C'l , there exists p m € i 1 (R w ) such that p m > 
a.e. and the projected measure p(t) — U#p,(t) satisfies 

p(t)({x})>0; 

c ) C' t ' is at most countable for each t e R; 
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d) if p m € L 1 (R JV ) is sucft t/iat p™ > o.e. on R N and p in = a.e. on J t _1 (R+) 
w/wte ||p*"||i,i(rjv) > 0, then 

p(t)({x}) > if and only if {p in ^ N ){F t ~\{x}) n ^({O}))) > 0. 

iYoo/. If a; € C t ", then Jgf ^(^({x}) n J t _1 ({°}) > so that in particular the 
intersection F 1 T l ({x}) n J t _1 ({0} is infinite, which implies that x g CtC\C' t , thereby 
proving a). 

Let R > 0; since F t is proper by Proposition 3.1, Ff (B(0,R)) — Kr is a 
compact subset of R w . Pick p m = 1k r ', obviously p m g L 1 (R Ar ) and p m > a.e.. 
Set p m = & N ®{l KR 5 uin ), let fj,(t) = $#p™ be the push-forward of p in under the 
Hamiltonian flow and let p(t) = H#p(t) be the projected measure. As explained 
above p(t) = F t #(l KR ^ N ) so that, whenever x g 5(0, fl) n C t ", 

= ^(F t - 1 ({x})) 

>^(F r 1 ({x})nj r 1 ({0}))>0. 

This immediately implies both b) and c). Indeed, b) is satisfied with p m = 1k r 
for each R > such that X g C" n S(0, J?). On the other hand, for each R > 0, all 
the points in C t " n -6(0, R) are atoms of the same Borel measure F t #(lK R £' N ) on 
R w , which implies c). 

As for d), one has again p(t) — F t #(p m Jz? N ), so that 

p(*)(M) = (p'^X^CM)) - (p'-^j^'tWlns) 

+ ( /) »^)(F ( - 1 (W)nJ t - 1 ({0})) 
+ (p i "if jr )(F t - 1 (W)nJ t - 1 (R;)) 

where £" is the Jz? ^-negligible set where F t is not diffcrentiable. The conclusion 
follows from the fact that ££ N '(E) = and p in = a.e. on J t _1 (R;), so that 

P(t)({x}) = { P in J? N ){F-\{x}) n Jr^W)) 

□ 

If p m g L 1 (R jv ) satisfies p ln > a.e. on K N and p in = a.e. on J t _1 (R+) 
while ||p m ||z,i(RN) > 0, one has (p m J£ N ) _L (J t ^f N ), so that the projected measure 
p(t) = H#p.(t) satisfies p(t) _L J£ N by statement b) in Theorem 5.1. By statement 
d) in Theorem 6.1, the measure p(t) is diffuse if and only if 

{p in J? N ){F t -\{x}) n J t -\{0}))) = for all x g C t . 

This may indeed happen, as shown by the following example. 

Example 5. Assume N = 1, and set H(x,{;) = ^£ 2 , which generates the free flow 
$t : (x, £) (a; + By regularity of the Lebesgue measure, there exists a 

compact set K C (0, 1) \ Q such that \ < ££ X (K) < 1. Let O = (0, 1) \ if; since ft 
is open in (0, 1) and contains (0, 1) n Q, it is a countably infinite union of disjoint 
nonempty open intervals: 

n = [J i n , so that jjf^n) = ^ ^(/n) • 

ngN riGN 
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In particular A := Jzf^fi) > since > for each n G N (indeed each I n is 

an open interval that contains at least one rational). Set t = 1 and 

(y, ify<o, 

^i(y) := < J V ln(z)dz, if ye (0,1), 
[ A + (y - 1) , if y > 1 . 

Define 

tT»(y) :=F(y)-y. 

The function Fi is Lipschitz continuous on R, being the indefinite integral of a 
bounded measurable function. Therefore U m is also Lipschitz continuous on R. 

The function F± is increasing on R — this being obvious on Rl and on (1, oo). 
Indeed, if < 2/1 < 2/2 < 1, the interval (3/1,2/2) contains at least one rational 
number r, so that 

F 1 (y 2 )-F 1 (y 1 )=J? 1 (nn(y 1 ,y 2 ))>0, 

since the open set fin (2/1,2/2) contains the intersection of the connected component 
of r in fi with (2/1,2/2), which is a nonempty open interval. Since F\{y) — > ±00 as 
2/ — > ±00, we conclude that F\ : R — >• R is one-to-one and onto. 

By the Lebesgue differentiation theorem (Theorem 7.11 in [22]), F\ is differen- 
tiable a.e. on R (in fact it is even of class C°° on H* U (1, 00)) and F[(y) = ln(y) 
for a.e. y G (0,1). Thus, denoting by E the ^f 1 - negligible set on which F\ is 
not differentiable, one has (-F 1 i')~ 1 ({0}) U E = K. Since F\ is one-to-one and onto, 
C x =Fx{K). 

Let p m = 1 K ; thus \\p in \\ L i = SS\K) > 0, so that (p™^ 1 ) _L (iFj^ 1 ). The 
measure p(i) = U#/j,(t) where /j(t) = $#/x m with p m := .S^ 1 <g> ^ n (x)5 uin{x) 
satisfies p(l) _L «S? N by statement b) in Theorem 5.1 

p(l)(R) = p(l)(Ci) = ||p in ||i 1 = if 1 (K)>0. 

On the other hand, since F\ is one-to-one and onto, for each x G R, one has 
#^i -1 ({ x }) = 1- I n particular, for each x G Ci = F\{K), one has 

FfHix}) C (Flr'dO}) so that #(F 1 - 1 ({x})n(F 1 )" 1 ({0})) - #F^{{x}) = 1. 
In particular, for each x G Ci, by statement d) of Theorem 6.1, one has 

( P ^ 1 )(Fr 1 ({x})n(F 1 ')- 1 ({0})) = 0, 

so that 

p(i)(W) = o. 

Hence p(l) _L ^f 1 and is diffuse, while p(l) ^ since p(l)(Ci) > 0. 

7. Application to the classical limit of quantum mechanics 

In this section, we apply the results obtained above to the classical limit of the 
Schrodinger equation. 
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7.1. The classical scaling. Consider the evolution Schrodinger equation 

ihd t ^ = -£ l ^ + V(x)^ 

for the wave function ^ of a point particle of mass m subject to the action of an 
external potential V = V(x) € R. 

Choosing "appropriate" units of time T and length L, we recast the Schrodinger 
equation in terms of dimensionless variables t := t/T and x := x/L. We define a 
rescaled wave function xp and a rescaled, dimensionless potential V by the formulas 

T 2 

4>(t, x) := Mt, x) and V(x) := —^V(x) . 

mL z 

In these dimensionless variables, the Schrodinger equation takes the form 
.„ 7 KT * mL 2 y 7 

The dimensionless number 2tt HT/mL 2 is the ratio of the Planck constant to mL 2 /T, 
that is (twice) the action of a classical particle of mass m moving at speed L/T on a 
distance L. If the scales of time T and length L have been chosen conveniently, L/T 
is the typical order of magnitude of the particle speed, and L is the typical length 
scale on which the particle motion is observed. The classical limit of quantum 
mechanics is defined by the scaling assumption 2t:H <C mL 2 /T — i.e. the typical 
action of the particle considered is large compared to h. Equivalently, mL/T is 
the order of magnitude of the particle momentum, so that 2irhT/mL is its dc 
Broglie wavelength; the scaling assumption 2nhT/mL <C L means that the de 
Broglie wavelength of the particle under consideration is small compared to the 
observation length scale L. 

Introducing the small, dimensionless parameter e = HT/mL 2 and dropping hats 
in the dimensionless variables as well as on the rescaled wave function and di- 
mensionless potential, we arrive at the following formulation for the Schrodinger 
equation in dimensionless variables 

(40) ied t ^ = -\e 2 & x ^ + V{x)i>. 

The WKB ansatz postulates that, at time t = 0, the wave function is of the form 
tp{t,x) = a in (x)e iSin W/ h , xeR". 
Consistently with the scaling argument above, we set 

a m (x) := a m (x) and S m (x) := TS m (x)/mL 2 
— since S m has the dimension of an action — so that 

(41) V>(0,z) =a in {x)e i§in W> e . 

Dropping hats in the initial data as well as in the Schrodinger equation, one ar- 
rives at the following Cauchy problem for the Schrodinger equation in dimensionless 
variables: 



(42) 



iedtipe = -^e 2 A x ip e + V(x)i(; e , x € R N , t G R , 
^ £ (0,x) = a in {x)e ism ( x)/e . 



The problem of the classical limit of the Schrodinger equation is to describe the 
wave function tp e in the limit as e — > + . 
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7.2. The WKB method. The traditional procedure for describing the classical 
limit of the Schrodingcr equation is the WKB method recalled below. First we 
recall some mathematical tools and elements of notation used in the presentation 
of that method. 

Assume that V <= C°°(R W ) satisfies 

(43) d a V e L°°(R N ) for each multi-index a E N N 
and 

(44) -> as Id -> +oo . 

fI 

Then the Hamiltonian 

#(*,0 :=\n 2 + V{x) 
satisfies (7) and therefore generates a global flow 

H N x R N 3 (x,0 i y {X t (x,t),Et(x,Q) E R N x R w 

by Lemma 2.f . 

Assume further that 



(45) 
with 

while 



sup / T v (x - y)V~{y)dy -> as 77 -> if JV > 2 

r„(*) = 



l [0M (\z\)\z\ 2 - N if TV > 3, 
l [0;?7] (M)ln(l/>|) if TV = 2 , 



/■x+l 

(46) sup / y-(y)dy<oo if N = 1 . 



Under assumptions (45)-(46), the operator —^e 2 A x + V has a self-adjoint exten- 
sion on L 2 (R N ) that is bounded from below. 

Under assumption (43), there exists a FIO that is a parametrix for the operator 

Ge (t) : = e'K^-iO 

see for instance Theorem 2.f in [17], whose main features are recalled below. 
Consider the action 

(47) S(t,x,Z):= [ (l\Z s (x,0\ 2 -V(X s (x,t)))ds 

Jo 

Given T > 0, we shall have to deal with the class of phase functions 

Lp = ip(t, x, y,r))eC of class C°° on [0, T)xK N xK N x R N ) 
satisfying the conditions 

' <f(t, X t (y, rj),y, 77) = S(t, y, 77) , 
D x ip(t,X t (y,rj),y,rj) = E t (y,r)) , 
iD 2 ip(t, x, y, 77) < is independent of x , 
k det( J D x ^(t,A t (7/, 77), 7/, 77)) ^ for each (t,y, V ) E [0,T) x K N x R N . 
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Pick x € C%°(R. N x R N ) and T > 0. Then, for any phase function (p satisfying 
(48) and any n > 0, there exists A n = A n (t,y,r),e) € C c °°([0,T] x R w x R")^] 
such that the FIO G n (t) with Schwartz kernel 

J<p(t,x,y,r,)/e_ df l 



G t , n (t,x,y) = J A n (t,y,r}, 



e e 



n-2N 



(2ne) N 

satisfies 

sup \\(G e (t) - G e , n (t)) X (x,-ied x )\\ c{L z {RN)) < C[V,T, X ]e 

0<t<T 

In this inequality the notation x( x , —ied x ) designates the pseudo-differential oper- 
ator defined by the formula 

X (x,-ied x ) ( f ) (x):= [[ e Kx- y Yvh x{x ^ {y) Wv 

Taking Theorem 2.1 in [17] for granted, one arrives at the following description 
of the classical limit of (42). It is stated without proof in Appendix 11 of [6] or as 
Theorem 5.1 in [4]. 

Let U in = VS in and let C be defined as in (15); let J\f(t,x) and y 3 (t,x) be 
defined as in Proposition 3.3 for each (t, x) e R x R N \ C. Let Jt(y) be defined as 
in (14). 

Proposition 7.1. Let a in e C™(R N ) and S in e C m+1 (R W ) with m > 6N + 4. 
For all e> and all (t, x) £ R + x K N \ C, set 

M(t> x ) in/ /+ \\ 

(49) * e (t,x)= J2 ^^e^^^, 
where 

Sj(t,x) := S^fafrx)) + S(t, Vj (t,x), \7S m (y 3 (t,x))) , j = 1, . . . ,A^,x) 

S(t,y,£) is given by (47) and v 3 (t,x) € Z /or a// (<,x) e R+ x R w \ C is constant 
on each connected component of R x R w where j < Af. 
Then the solution ip e of the Cauchy problem (42) satisfies 

(50) ip e (t,x) = V e (t,x) + R\ (t,x)+R 2 e (t,x) 
for all T > 0, where 

SU P \\RI\\l 2 (b(o.r)) = 0(e) for all R > 

0<t<T 

and 

sup \R 2 (t,x)\ = 0(e) for each compact K C R+ x R N \ C 

(t,x)eK 

as e -» + . 

A self-contained proof of Proposition 7.1 based on the parametrix construction 
of [17] is given in Appendix A . 

The integer u 3 (t,x) in (49) is defined precisely in the proof (see formula 66 
below). How this integer is related to the Morse index of the path 

[0,t] 3 s^X s (y 3 (t,x),VS m (y 3 (t,x))), 

is not needed in the sequel and will not be discussed here. We refer to [4] for more 
details concerning this point. 
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7.3. The Wigner transform. Alternately, the asymptotic limit of ip e can also be 
investigated with the help of Wigner 's transform [26, 18]. 

For each ^ e L 2 (R, N ), one defines the Wigner transform of ^ at scale e by the 
formula 

W e [*]{x,Z):= [ e-^(x + \ey)W^hy)7TTN- 

J R N {ITT) 

Let us define the e-Fourier transform T t as follows: 



(R») 

Given b = b(x,t;) in the Schwartz space S(R N x R N ), we recall that the e- 
pseudodifferential operator B e with Weyl symbol b is the integral operator with 
integral kernel 



X + V A M-(x-y)/e dtj 



Straightforward computations lead to the following identities: for each * e 
S(R N ) 

(51) f W e [*](x,£)dS = l*(*)| 2 , / W e [*](x, £)dx - \^(0\ 2 ■ 

Likewise, for each e L 2 (R N ) and each e-pscudodiffcrential £? e operator with 
Weyl symbol b(x,t;) in the Schwarts space ^(R^ x R"), one has 

(52) f W e [^}(x,0b(x,0dxd^ = (^,B e ^) L 2 {nN) . 

Lemma 7.2. Let a in E L 2 (R jv ) and S in e ^(R"). Ifcen 

VF £ [a m e lS *"/ e ](a;,0 -> (a in ) 2 (a;)<K£ - V5 m (x)) in ^(R^ x R w ) 

as e — > 0. 

See Example III. 5 in [18]. 

Proposition 7.3 ([18]). Assume that V e C^^R") satires (45) -(46) -(43) -(44). 
Let 

for each e > and eac/i t € R, and /or a.e. (x, £) G R w x R w . 

T/ien i/ie family /U e (t) converges in S'(R N x R w ) uniformly in \t\ < T for each 
T > to i/ie unique solution of the Liouville equation 

jd tt i + t- v x/ u - v x y (x) ■ v £ m = o , x, e e R w , t e R , 

\MO,^,0 :=a <n (x) 2 5 vs4 „(x)(0- 



(53) 



in particular, \i e C;,(R, w — M(R N x R^)) and is a positive Radon measure 
for all t G R satisfying 



J Jr»xr» 



/or eac/i t e R. 
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This proposition is stated as Theorem IV. 1 in [18], to which we refer for a 
complete proof. 

Therefore, the structure of the measure fj,(t) is given by Theorem 5.1 with 

p m = { a m f and U m = V S m . 

Likewise, Theorem 6.1 provides a precise description of the exceptional sets associ- 
ated to fl(t). 

7.4. WKB vs. monokinetic measures: the case of rough phases. The WKB 

asymptotic solution (49) obviously contains more information than the regular part 
/U a (i) of the Wigner measure fi(t) that is the solution of the Liouville equation (53) 
— for the definition of /i a (t) and its structure, see Theorem 5.1. 

For instance, the (finitely many) phase functions Sj(t, x) appearing in (49) define 
velocity fields Uj(t,x) :— V x Sj(t, x) as in the formula giving fx a (t) in Theorem 5.1. 
On the contrary, the nonnegative integers Vj{t,x) do not appear in the formula for 

One can go a little further, and compute the Wigner measure associated to the 
WKB asymptotic solution (49), under the assumptions used in Proposition 7.1. 
Let t e R; for each x e R w \ C t and j € {1, . . . ,Af(t, x)} 7 one has 

V x Sj(t,x) = DS in ( yj (t,x))D xyj (t,x) + D y S(t, yj (t,x),DS in ( yj (t,x)))D xyj (t,x) 

+D v S(t, Vj (t, x),DS m ( Vj (t, x)))D 2 S m ( yj (t, x)D xVj (t, x) 

= E t ( yj (t,x),DS in ( yj (t,x))) ■ (D y X t ( yj (t,x),DS in ( yj (t,x)))D xyj (t,x) 

+D n X t ( yj (t, x),DS m (y,(t, x)))D 2 S m (y,(t, x))D xVj (t, x)) , 

in view of formulas (3.1-2) in [17] recalled above (see (61)). By definition 

(54) X t ( yi (t,x),DS m ( yj (t,x))) =x, j = l,...,M(t,x), 
so that 

D y X t { yi {t, x), DS m ( Vj (t, x)))D x yj(t, x) 
+D v X t (y 3 (t,x),DS m (y 3 (t,x)))D 2 S m (y 3 (t,x))D xyj (t,x) = I 
which implies in turn that 

(55) W x Sj(t,x) = ~ t ( y] (t,x),DS m ( y] (t,x))), j = l,...,Af(t,x). 
Because of (54) and of the condition 

yj(t, x) ^ yk (t, x) if 1 < j < k < Af(t, x) , 

one has 

V x Sj (t,x)= St ( Vj (t,x), DS m ( Vj (t,x))) 

^Z t (y k (t,x),DS m ( yk (t,x)))=X7 x S k (t,x) if 1 < j <k<M(t,x), 

by uniqueness of the solution of the Cauchy problem for Hamilton's equations. 
Applying Proposition 1.5 in [10] implies that 

W e [* e (t, •)] m in V'((R N \ C t ) x K N ) 

as e — > 0, where 

^ V ^' Jt(Vj(t,x)) v * bAt ' x) 
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The expression in the right hand side of the equality above coincides with the 
regular part fi a (t) of the Wigner measure fi(t) in Theorem 5.1. 

Obviously Theorem 5.1 (bearing on Wigner measures) holds under assumptions 
on the phase of the WKB type initial condition much weaker than in Proposition 
7.1 (bearing on solutions of the Schrodinger equation). How much information 
on solutions of the Schrodinger equation can be extracted from Theorem 5.1 is 
therefore a natural question. Since the Wigner measure is the limit of the Wigner 
transform as e — > and the Wigner transform is quadratic in the wave function, it 
is natural to expect that Theorem 5.1 carries information on the vanishing e limit 
of appropriate quadratic expressions in the solution of the Schrodinger equation 
(42) with rough initial data. 

Using (51) and (52), we arrive at the following statement. 

Proposition 7.4. Let ip e the solution of the Schrodinger equation (42) where a m G 
L 2 (R N ) satisfies ||a m || L 2 (RN) = 1 and S m G C^R") is such that U m := V5 m 
satisfies (11) and the regularity condition (18). Assume that the potential V G 
C 1 ' 1 ^) satisfies (45)- (46)- (43)- (44). Then 

a) for each e-pseudodifferential operator B e of Weyl symbol b = b(x,£) in the 
Schwartz class ^(R^ x R^) and for each t G R, one has 

lim(ip e ,B e ip e ) L 2 {R N ) = / b(x,£)fj,(t,dxd£) 

e^O Jr»xR» 

where /i(t) is the Wigner measure described in Theorem 5.1; 

b) for each t G R and each x <= ^(R^) satisfying 

X{x) — for all x G C t , 

one has 

[ x{x)\Mt,x)\Hx^ [ X (x) ]T |ffl ' n|2lpt (y)dx 

yeF t 1 ({ X }) 

as e — > 0; 

c) for each t E R and each x G C^R*) satisfying 

X (Zt(y,VS m (y)) = for all y G F^{C t ) , 

one has 

[ X (0\^Mt,0\ 2 d^ f E x(Zt{y^S m {y))) ]am] \ 1Pt {y)dx 
Jn N Jr« J~r., J t 

as e — > 0. (Recall that C t is defined in (27) — see also (22) and (19) — while P t , 
Jt, S t and F t are as in Theorem 5.1.) 

Compared with statement a), statements b) and c) involve a tightness issue, 
i.e. the fact that W e [tp e (t, •)] is small at infinity uniformly in e G (0,1). This is 
discussed in [18] after the statement of Theorem IV. 1; see also the end of the proof 
of Theorem 2.1 on p. 337 in [10]. We give the proof of Proposition 7.4 below for 
the sake of being complete. 
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Proof. Statement a) follows from formula (52) and Proposition 7.3. 
Now for statements b) and c). Denote 

:= J^WTi^ 2 ^ > G *( X >® ■■= 9e(x)9e(0 ■ 

Defining 

** <«°(x) = ^ e yN/i e -\x- Xo \ 2 /2e e iio-x/e _ 

one has 

Along with the Wigner transform W e [i[) e (t, ■)], we consider the Husimi transform 

Applying (52) with b — G e and ^ = t/> € (i, •), one finds that 

e (2^e) wl e M e 1 
so that, for each £o>£o € R w and each e > 

W e (t, ■)](*«,, 6>) - 7^pvK^ ' l0 l^^-))| 2 > 0. 
Elementary computations based on this formula show that 

/ W t [ip e (t, ■)](*(>, - (.9, * |Ve(V)| 2 )M 

(56) Y _ 

/ W e {Mt, ■Mxoi&dx = (g e * \F e 4>e{t, -)| 2 )(Co) 

so that 

// W e [ip e (t,-)](x ,£o) = // 5e(a; - 2o)|-0 e (i,2O| 2 d£cteo 

JJr»xr» J Jn N n N 

= IIV>e(V)||l,2(R*) = ll^rilL2(R«) = ||a m ||L2(R«) = 1. 

On the other hand 

// n(t,dxd£)= [[ f i(0,dxdO = \\a m \\ L 2 {RK) = l 

JJR«xR« JVxR" 

since /z(f) is the push-forward of the probability measure /u(0) under the Hamilton- 
ianflowof ±|£| 2 + V(x). 

Since W e [ip e {t,-)] > and 

(57) ff W e [Mt,-)]{x,0x{x,0dxdS^ ft X (x,0^(t,dxdO 
for each X € C C (R JV x R*) while 

W^t, -)](x,0dxdS = l=H »(t,dxdO , 

we conclude that (57) holds for each \ <E Cb(R N x R w ) (see for instance Theorem 
6.8 in chapter II of [19]). 
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On the other hand, for each ^eCj (R ) 

/ X (x)(\Mt,x)\ 2 -\Mt,-)\ 2 *9 e (x))dx 
Jr n 

< \x{x) - X*x gAx)\\ipe(t,x)\ 2 dx 
<\^||V X ||l~ I \y\gi(y)dy ^ 



and likewise 



0. 



We conclude that, for each \ g Cl(R N ) 



(58) 



/ X (x)\Mt,x)\ 2 dx^ // X {x)n(t,dxdO 

JR» iJR s xR» 

/ X (0\^Mt,0\ 2 d^ II x(ZMt,dxdt) 

JR» JiR»xR» 



as e — » 0. On the other hand, 



\ll>e(t, x)\ 2 dx 



|^ £ (t,0l 2 « 



R"xR N 



/j,(t, dxd£) 



so that (58) holds for each x S CjfR^) by a standard density argument. 

Statement b) follows from the first convergence statement in (58): by Theorem 
5.1 e), 

// x(x)v(t,dxd£) = x(x)p(t,dx)= X (x)pa(t,dx) 

JJr«xr» Jr« Jr." 

if the test function x vanishes on C t . One concludes with the formula for p a in 
Theorem 5.1 d) with p m = |a m | 2 , the restriction to H N \C' t being unessential since 
& N {C' t )=Q. 

On the other hand, if x(S t (y, VS in (j/))) = for all y g F^{C t ), then 
/ x{t)n{t,dxd£) = I x(0 1 n N \c t (x)p(t,dxdO 

JR N xR N iR«xR N 

by Theorem 5.1 c). By Theorem 5.1 b) and e), since n(t) is a positive measure, its 
singular part n s {t) is carried by C t x R N , and therefore 



L 



R"xR» 



x(0 1 B. N \c t (x)»(t,dxd4) 



R"xR» 



xK)1r«\c 4 (x)n a (t,dxd4) 



r \ n tn\A-t 

/ Yl x(Zt(yyS m (y)) l -^-^(y)dx. 



This proves statement c). 



□ 
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7.5. Conclusions and perspectives. Although Proposition 7.1 can be improved 
by adding extra assumptions on a m and S m , in full generality the approximation 
stated there holds in R x R N \ C, and only after localization in compact subsets of 
R x R N \ C. In particular, the approximation of the wave function ijj e by the WKB 
ansatz \I/ e is not uniform as (t, x) approaches C. (For more on this well known fact 
in a different, yet related context, we refer the reader to the last paragraph in §55 
and on the problem in §59 in [15].) 

On the contrary, the result obtained in Theorem 5.1 is global, holds in the sense 
of measures and does not involve the caustic C as explained in the previous sections. 

Even in the case of rough phase functions, statements b) and c) of Proposition 7.4 
show that our approach provides information on the solution i/j e of the Schrodinger 
equation(specifically, on \tp e \ and |J r e '0e|) consistent with a WKB ansatz, in a sit- 
uation where the approximation (50) by a WKB asymptotic solution of the form 
(49) is not justified at the time of this writing. 



Let x(x,t) = Xi(x)X2(0 with xi,X2 & C™(R N ), satisfying 

1b(o,r)(x) < Xi(x) < 1b(o,_r+i) and l B (o,Q)(£) < X2(0 < lfi(o,Q+i) 
for all x, £ <£ R N , where R > and Q is to be chosen later. Pick n > 27V; then 



for all t e [0, T], where Cq.t = C[V, T, X }- 

Now, x(x,—ied x )tp l e n = xi(x)X2((-i£d x )ipl n and since supp(o* n ) C B(0,R) 



Appendix A. Proof of Proposition 7.1 



Mt, •) - G e , n (t) X (x, -ied x )^T - Ge(t)(l - X (x, -ieflxM" 

+ (G e (t)-G e , n (t)) X (x,-ied x )il>r. 



Since G e (t) is a unitary group on L 2 (Il N ) 



\\Mt,-) - G € , n (t) X (x,-ied x )ti n \\ LHRN) < \\(l-x(x,-ied x W £ n \\ LHIlN) 

+\\(GS)-G e , n (t))x(x,-ied x )4 n \\ L 2 {RN) 
< \\(l- X (x,-ied x ))4 n \\ L 2 {RN) 




+Ct,q£ 



||(l-X(^-^))Vni^(R-) = IIXi(l-X2HeO x M n |l^(R-) 

<||(l-X2He5 x M n H^(R-) 
= (2^-^11(1 - X2(eO)#1^(R") 
< (2^)- Ar ||l [Q/e , oo) (|£|)^F|| L2(R « ) 



Since 




we conclude from estimate (7.7.1') in [12] that 




.m 
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provided that supp(a m ) C B(0,R) and |C| > 1 + ||V5 J 



|L°°(S(0,K) 



. Therefore 



(59) 



\\1> e {t,-)-G e , n {t) X {x,-ied x )il>?\\ 



L 2 (R«) 



< (J _ 1 1 I _ ... ^n—2N 



J T,Q\\a ||L2( R N)£ 



+c\\(i + \c\r m \\ L ^M 



for all to > N/2. 

Next we analyze the term 
(60) 



A n (t,y,r],e)a m {z)x2(()e 



i(<p(t,x,y,r ] )+(-(y-z)+S* n (z))/. 



dzd(dydri 
(2Tre) 2N 



(61) 



with the stationary phase method. 
Choose ip of the form 

'Pit, x, y, rj) = Sit, y, r)) + (x- X t (y, i])) ■ E t (y, r))+iB : (x - X t (y, r/))® 2 

where the matrix B = B T > is constant (see formula (2.7) [17] and the following 
Remark 2.1). Critical points of the phase in the oscillating integral (60) are defined 
by the system of equations 

' -( + DS in iz)=0, 
y- z = 0, 

d y S{t,x,y) - E t (j/,T7) • D v X t {y,rf) + {x - X t {y,rf)) ■ D y E t (y,ri) 
- iB : (x - X t {y, rf)) ® D y X t {y, rf) + Q = Q, 
d v Sit,x,y) - H t (j/,T7) • D v X t (y,rj) + (x - X t (y,r])) ■ D ri E t (y,r]) 
-iB:(x- X t (y, jy)) ® D v X t (y, r,) + ( = . 
At this point, we recall formulas (3.1-2) from [17] 

d y S(t, x, y) = E t (y, rf) • D y X t (y, rf) - r] , 
dr,S{t,x,y) = Z t (y,rj) ■ D v X t (y,r]) , 
together with the following definitions 

Y(t, V, V) : = D y-tiy, rf) - iBD y X t (y, rf) , 
Z(t,y,rj) := D,S t (j/,jy) - iBD v X t (y,r]) . 
Thus the critical points of the phase in (60) are given by 
' ( = DS in iz), 
V = z, 

(x-X t (y, V )) T Y(t,y, V ) + ( = i 1 , 

[ix-X t (y,r])) T Z(t,y,r)) = 0. 

Since the matrix Z is invertible by Lemma 4.1 of [17], we conclude that the system 
of equations above is equivalent to 

C = DS m {z) , 

y = z, 
C = v, 

x = X t (y,rj) . 
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j = l,...,Af(t,x). 



In other words, 

F t (z) = x, y = z, C = V = DS m (z). 
Assuming that (t, x) £ C, we apply Proposition 3.3 and conclude that the set of 
critical points of the phase in (60) is of the form 

j y = z = yj(t,x) , 

\C = V = DS in ( yj (t,x)), 

At this point, we apply the stationary phase method (Theorem 7.7.5 in [12]). First 
we need to compute the Hessian of the phase in (60) at its critical points. One finds 



Hfax) := 



V = Vj (*,*) 



/ D 2 S in 


-I 





° \ 


-I 





+1 








+1 


-YD y X t 


-YD v X t - I 


V o 





-ZD y X t 


-ZD v X t J 



and it remains to compute det(Hj(t,x)). Adding the first row of Hj(t,x) to the 
third row, one finds that 



det Hj(t, x) 



V = Vj (t,x) 
,= DS'"( s ,( t|I )) 



V=Vj (*.*) 
>|=DS'»( Sj (t,.)) 



D 2gin 


-I 








-I 





+/ 





D 2 S" 





-YD y X t 


-YD, n X t - I 








-ZDyX t 


-ZD v X t 




-I 


+1 







D 2gin 


-YDyX t 


-YD, n X t - I 







-ZDyX t 


-ZD v X t 



= (-1) 



N 



-I 
D 2 S in 







-YDyXt + D 2 S m 

-ZDyX t 





-YD v X t - I 
-ZD v X t 



V=Vj (t,a=) 



where the last equality follows from adding the first column in the right hand side 
of the second equality to the second column. Eventually, one finds that 

-YDyXt + D 2 S m —YD v X t - I 
— ZDyXt —ZD. q Xt 



det Hj(t, x) 



which is computed as follows. First 
/ -YZ- 1 



/ 

so that 



-YDyXt 



D 2 S in —YD v X t 



-ZDyXt 



-ZD v X t 



det Hj(t, x) 



D 2 S in _j 

—ZDyXt —ZDfjXt 



y=Vj 

l)=DS'»( S ,(i,l)) 



D 2 S in _j 

—ZDyX t —ZD. q Xt 



V = Vj (t,x) 
■,=DS»( Sj (t,i)) 



On the other hand 

D 2 S in 



ZDyXt —ZD v Xt 



det(ZD y X t + ZD n X t D 2 S in ) 



det(Z) det(D y X t + D v X t D 2 S in ) 



by the following elementary lemma (that is a variant of the Schur complement 
formula in a special case: see for instance Proposition 3.9 on pp. 40-41 in [23]). 
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Lemma A. 1. Let A, B, C, D e M N (C). If AB = BA, one has 

A B 
C D 



= det{DA - CB) . 



(62) 



Therefore 

det Hj(t, x) = dct(Z) dct{D y X t + D v X t D 2 S in ) 



y=y j (t,x),r]=DS in (y j (t,x)) 



det(Z( Vj (t, x), DS m ( Vj (t, x))) det(DF t ( Vj (t, x))) 



where F t is defined in (13) with U in = DS in . 

Pick a nonempty closed ball BcRx ~R N \ C, let As = Af(t, x) for all (t, x) E B, 
and let 

Kj = {(y J (t,x),VS m (y 1 (t 1 x))) \ (t,x) e B} , j = l,...,Af B . 

Assuming that B is of small enough radius, Kj (~l Kk = for j ^ k £ {1, . . . , A"b}- 
Let Kj e C^°(R 27V ) for all j = 1, . . . , As, such that 

{Kj > and Kj\ K = 1 , j = 1, . . . , Mb , 
while KjKk — for j ^ k e {1, . . . , As} , 

Applying Theorem 7.7.1 in [12] shows that 
(63) 

'^(t, 2/ , ?? , e ) a »(z) X2 (C) e ^( t --'" + ^-^ +SS ' , ( z »/ e ^|^ 



sup 

{t,x)eB 



Hill Mt,y,V,e)a in (z) X 2(0^(y,v)^(z,0 



i( ip (t,x, v , v )+c-(y-z)+s>"(z))/e dz d(dydri 

(27re) 2W 



= 0(6) 



as e — > 0. 
Next we set 



Ij(t,x,e) := jjjj A n (t,y,r 1 ,e)a m (z)x2(0K J (y,r 1 )K J (z,C) 

yr t( V (t,x,y,v)+C-(v-z)+S'"(z))/e dzd C d y d V 

for j = 1,. . . , As- By Theorem 7.7.5 in [12], we conclude that 
(64) 

sup \Ij (t, x, e) - A (t, yj (t, x) , \7S m (y, (t, x)), 0)a m ( Vj (t, x)) X 2 (VS m ( Vj (t, x))) 
(t,x)eB 

xe <( v (( W (W,vs , ^ w (t,i)))+^(y J M))/ e ^ ff .^ J jj-i/a| = 

as e — > 0. Our choice of %2 and 95 implies that S m {yj{t 1 x))) = 1 and 
^(t,x,yj(t,x),VS in (yj(t,x))) = S(t )!/j (i ) a ! ) ) VS in (i/ j (i ) a!))) 

so that 

fo>(t, x, W (t, x), VS m (y 3 (t, x))) + S in ( yj (t, x))) = Sj(t, x) . 
By formula (2.13) in [17] 

(65) A Q (t,y ] (t,x),VS m (y J (t,x)),0) = ^ Det(Z( Vj (t,x), DS™ (y 3 (t,x))), 
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where the notation con ^fz designates the analytic continuation of the square-root 
along the path t i->- Det(Z(yj(t,x),DS m (yj(t,x)). This analytic continuation is 
uniquelcy defined since Det(Z(yj(t, x), DS m (yj(t, x))) ^ for all t: see for instance 
section 1.3 in chapter 8 of [1]. According to (65), (62), Lemma 5.1 and formula 
(5.15) in [17], we have 

A (t, Vj (t, x), VS m {y 3 (t,x)),0) (dct Hj (t, x))- 1 ' 2 

(66) 



|det(DF*( % (M)))r 1/2 e 1 ™^/ 2 = J( % (M)r l7 V 



■ Vj (t,x)/2 



where fj(t,x) is an integer. 

Putting together (59)-(63)-(64)-(66) concludes the proof of Proposition 7.1. 

Proof of Lemma A.l. If A is nonsingular and AB = BA, one has 



= (-1) 



N 



A- 1 




B 

-A 



= (-1) 



N 



I 

CA- 1 





CB-DA 



= (-1) N det{CB - DA) = det(DA - CB) . 

Since both sides of the identity above are continuous functions of A and the set 
of nonsingular matrices GLjv(C) is dense in Mjv(C), this identity holds for all 
A e M„(C) such that AB = BA. □ 
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